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ABSTRACT 

N-body simulations show that when infall reorientates the outer parts of a galactic halo 
by several degrees per Gyr, a self-gravitating disk that is embedded in the halo develops 
an integral-sign warp that is comparable in amplitude to observed warps. Studies of 
angular-momentum acquisition suggest that the required rate of halo reorientation is 
realistic for galaxies like the Milky Way. 
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1 INTRODUCTION 

Galactic warps are common but remain mysterious. Bosma 
(1978) showed that of order one half of all galactic HI disks 
are measurably warped, as is the disk of the Milky Way. 
Briggs (1990) studied the warps of a sample of 12 galaxies 
in some detail and inferred several general laws that govern 
the phenomenology of warps. 

Theoretical studies of warps have made substantially 
less progress, so that after four decades of work on this 
problem, there is no concensus as to how warps work and 
why they are so common. An explanation that has enjoyed 
some success since it was first mooted (Toomre 1983; Dekel 
& Shlosman 1983) is that warps persist indefinitely because 
they are normal modes of oscillation of the disk in the poten- 
tial of a flattened massive halo. Recently it has emerged that 
this proposal fails completely when proper account is taken 
of the internal dynamics of the halo (Nelson & Tremaine 
1995; Binney, Jiang & Dutta 1998 - hereafter Paper I). Con- 
sequently, there has been renewed interest in other possible 
explanations. Ostriker & Binney (1989) pointed out that 
the net angular-momentum vector of a typical galactic halo 
must significantly shift its direction every Hubble time or so, 
and that the resulting torques are liable to warp any galac- 
tic disk. In particular they investigated the possibility that 
such torques caused the misalignment of the Milky Way's 
inner and middle disks. 

Ostriker & Binney did not attempt to reproduce the 
Galaxy's warp in the same way because of the undefined 
nature of the dynamics of the dark halo in the case that 
the halo is disturbed by accretion. Garcia-Ruiz, Kuijken & 
Dubinski (1998) briefly report a first attack on this prob- 
lem, in which the halo is rotated rigidly. Here we describe 
experiments in which a fully dynamical halo of 10 5 particles 
reorientates its angular-momentum vector as a result of ac- 
creting fresh material. A disk comprising 100 massive rings 
is embedded in the halo. We show that this disk develops an 
integral-sign warp very similar to those actually observed. 

2 THE COMPUTER CODE 

Our computer code and model galaxy are described in de- 
tail in Paper I, so here we only summarize them and indi- 



cate differences between the present simulations and those 
described in Paper I. The halo is represented by a standard 
A r -body code in which the potential is found by expanding 
the density on a polar grid in spherical harmonics. The disk 
consists of massive rings that couple to each other by di- 
rect summation of the vertical forces between rings, under 
the assumption that the maximum excursion of any ring 
from the xy plane is small compared to the ring's radius. 
The torques on rings from the halo are calculated from the 
spherical-harmonic expansion for the halo's potential. Con- 
versely, the rings exert forces on the halo particles because 
the rings contribute to the latter expansion. Thus, com- 
plete symmetry is maintained in each type of interaction 
(ring-ring, ring-particle and particle-particle) while ensur- 
ing that the disk is a tightly coupled entity. The sum over 
spherical harmonics is carried up to and including I = 8 and 
then truncated. Paper I describes numerical checks on the 
accuracy of the code. 



3 THE GALAXY MODEL 

The surface density of the disk, which is represented by 100 
rings, is exponential with scale-length R& out to Rt = 3.5Rd- 
At R > R t the surface density is smoothly tapered to zero 
at Rq — 4Rd- At the start of the simulation the disk lies 
flat within the equatorial plane of the halo. 

The halo is represented by 10 5 particles and has a total 
mass ten times that of the disk; consequently, the average 
mass of a ring is 100 times the mass of a halo particle. The 
particles are distributed such that the circular-speed of the 
combined disk and halo is very nearly flat out to 5Rd- Fur- 
ther out the circular-speed curve gradually drops, becoming 
Keplerian at R > 14Rd- The halo is flattened along the 
symmetry axis of the disk, with initial axis ratio ~ 0.77. 
In one simulation the halo has no net rotation. In another 
the halo is caused to rotate in the same sense as the disk 
by reversing the angular momentum, L z , about the disk's 
spin axis of one fifth of the halo particles: this was done by 
picking halo particles at random and reversing L z if this was 
opposite in sign to L z for the disk. 
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Figure 1. In each panel the curve shows, as a function of radius, the inclination 8^ of the disk to its original symmetry axis while the 
points show the inclination of the halo to the same axis. Filled symbols indicate that the halo is inclined in the same sense as the disk, 
and open symbols indicate the opposite sense of inclination. Beyond the edge of the disk, the halo's inclination is marked by triangles. 
Between successive panels a time has elapsed which is equal to ~ 2.3 times the time, t re f, required to travel the halo's half-mass radius 
at the rms speed of a halo particle. In this simulation the halo has no net rotation. 



4 MODELLING INFALL 

We model infall of material by steadily injecting particles in 
a toroidal region. Specifically, in an appropriately oriented 
coordinate system, each added particle has coordinates 

x = (a + s cos 9) cos 4> 

y = (a + s cos 9) sin <f> (1) 
z = ssin#, 

where a — 8.9i?d and 9 and 4> are angles that are uniformly 
distributed in (0, 2ir). The variable s in equation (1) is a 
random variable chosen such that the surface probability 
density is 

P(s) = ±e-° 2 ^ 2 (2) 

with 6 = 0.11i?d- Thus, the injection torus is character- 
ized by radii a and 6 with a/b = 80 and has an approxi- 
mately Gaussian density profile when sliced perpendicular 
to its larger circle. The symmetry axis of the torus is tilted 
by 15° with respect to the disk's initial spin axis, so that the 
torus passes above the x axis of the disk at x > 0. The ini- 
tial velocities of injected particles are directed tangentially 
around the torus, in the sense of disk rotation, and are equal 
in magnitude to the local circular speed at the moment of 
injection. 



A useful characteristic time is the time required to 
travel a distance equal to the halo's half-mass radius (~ 
3.2i? d ) at the rms speed of a halo particle. If we call this 
time t Te t, then a fresh batch of particles was injected each 
2.3t re f such that the number of particles that had been in- 
jected up to time t is given by 

N ini (t) = 15 000 [l - cos(7rf/t max )] , (3) 

where t max = 25t rc f is the duration of the simulation. As 
equation (3) implies, in all 30 000 particles, or 30% of the 
original halo mass, were injected. 

5 DIAGNOSTICS 

We wish to describe the orientation of each of the disk's 
rings and to relate this to the shape of the halo's mass dis- 
tribution at a similar radius. To this end, for each ring, 
we determine the angle 9i between the ring's axis and the 
original symmetry axis of the disk, as well as the azimuthal 
angle <f>i of the ring's lowest point. 

We determine the shape of the halo as a function of 
radius r from the coefficients Ai m (r) and Bi m (r) used to 
expand the halo's density in spherical harmonics [equation 
(16) of Paper I]. As is well known, for given I the complex 
coefficients Xi m = Ai m — iBi m transform under rotations 
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Figure 2. The same as Fig. 1 except that the angles marked are the values of <j>d (curves) and <j>h (points). 



of the coordinate system that defines the polar angles, as 
the spin-/ irreducible representation of the rotation group, 
0(3). Hence, once the Xi m have been determined in any 
given frame, their values X' lm in any rotated frame may be 
obtained by multiplying by the spin-/ matrix for the rela- 
tive orientation of the two frames. The coefficient X' 20 (r) 
will be largest when the z axis of the rotated system most 
nearly coincides with the shortest principal axis of the halo 
at radius r, and we determine the direction of this axis by 
maximizing IX20I 2 with respect to the Euler angles 9 and <f> 
upon which it depends. Since this procedure does not dis- 
tinguish between a direction and its opposite, it does not 
distinguish between (9, <f)) and (9' , 4>') = (tt — 9, (f> + tt). We 
resolve this ambiguity by requiring 9 to be smaller than ir/2. 

6 RESULTS 

Figs 1 and 2 compare, as functions of r, the polar angles 9d 
and (j>d that specify the direction of each ring's symmetry 
axis (curves) with the angles 9 h and 4> b that specify the 
direction of the shortest principal axis of the halo (points). 
For ease of interpretation, the halo's angles are plotted as full 
points when \<j>d — 4>h\ < 90°, and as open points otherwise. 
Hence, full points indicate that the halo and disk are tilted 
in similar directions, while open points indicate that the two 
components are tilted in opposing directions. 

In the top left panel of Fig. 1, which describes the initial 
condition, the disk's curve is invisible because it runs along 



the x-axis and the halo points give an idea of the statistical 
error in the determination of 9^. In the second panel the 
disk can be seen to have developed a small corrugation. The 
equal numbers of open and filled squares that describe the 
halo in this panel indicate that the disk and the halo are 
not systematically aligned at this time, probably because 
the disk's small corrugation merely reflects statistical noise 
in the simulation. By the panel labeled 3, the outer halo 
has clearly developed a systematic tilt of order 3°. By the 
time of panel 5, the preponderance of filled halo points in 
the region of the disk suggests that the halo's tilt has been 
to some extent communicated to the disk. In the next few 
panels it becomes steadily clearer that the tilt is propagating 
from the outer halo inwards, and causing the inner halo and 
disk to tilt very much as a unit. By the time of panel 9 both 
9h and 9d increase essentially linearly with radius from ~ 2° 
at r < Ed to ~ 5.5° at r ~ 4_R d - The increase of 9 d with 7? 
implies that the disk has developed an integral-sign warp, 
with the outermost ring moving ~ 0.3i?d out of the plane of 
the innermost ring. For the case of the Milky Way, for which 
i? d ~ 3kpc, this would correspond to a warp amplitude of 
order 1 kpc at 14 kpc, very similar to that actually measured 
(Burton 1992). 

Fig. 2 compares the azimuth of the ascending node of 
the disk, <j)d, with the corresponding azimuth, 4>h for each 
spherical shell of the halo. Initially, <f>d is arbitrarily set 
to 90° and statistical fluctuations cause <^> h to be widely 
scattered. By the time of panel 5, the scatter in 4>h has 
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been reduced to ~ 40°, and this scatter continues to de- 
crease as further infall causes the halo to become more and 
more inclined to the z-axis of the coordinate system. The 
disk is closely aligned with the halo from panel 5 on. The 
mean value of <f>d shows a statistically marginal tendency 
to decrease from the times of panel 8 to that of panel 11, 
consistent with the disk precessing retrogradely in the grav- 
itational field of the outer halo. It is not clear, however, 
that this precession, if statistically significant, is other than 
a transient phenomenon. 

The results for the case of a slowly rotating halo are 
essentially identical to those for the non-rotating case that 
are shown in Figs 1 and 2. 

7 CONCLUSIONS 

We have shown that a realistic warp develops in a disk em- 
bedded in a live halo when the halo accretes material whose 
angular-momentum vector is inclined slightly with respect 
to the original symmetry axis of the galaxy. Is this mecha- 
nism the long-sought solution to the warping problem? 

The close alignment of the disk with the local halo that 
is evident in both the simulations presented here and in Pa- 
per I indicates that a warp is essentially a halo phenomenon: 
the disk serves largely as a handy tracer of the halo's dy- 
namics. The essential feature of our halo simulation is the 
realignment of the approximate symmetry axis of the halo 
at r ~ 10i?d by ~ 7° over the duration, i max , of the simula- 
tion. If we scale the simulations to parameters characteristic 
of the Milky Way (R d ~ 3kpc, v c ~ 220 kms -1 ), we find 
i max ~ 0.9 Gyr. Hence, we would seem to have a convincing 
explanation of warp if the symmetry axis of the outer halo 
can shift by ~ 7° in ~ 0.9 Gyr. 

Several authors have studied how tidal interactions en- 
dow protogalactic regions with angular momentum at early 
times (e.g. Heavens & Peacock 1988) and these studies are 
generally in good agreement with numerical simulations 
(Barnes & Efstathiou 1987). Ryden (1988) and Quinn & 
Binney (1992) investigated the rate at which the direction of 
infalling angular momentum should slew by studying the an- 
gular momenta of individual spherical shells of protogalactic 
material. Quinn & Binney found that the angular momenta 
acquired by shells that differ in radius by a factor 2 have a 
clear tendency to be antiparallel. Moreover, the angular mo- 
mentum per unit mass of a shell rises strongly with radius, 
with the consequence that the net angular momentum of a 
galaxy tends to be aligned with the angular momentum of 
the most recently accreted shell. These two results together 
imply that the net spin axis of a halo tends to slew through 
more than 90° in the time required for the radius of the 
currently accreting shell to double. This time depends on 
the cosmology and the initial density profile (e.g. Fillmore 
& Goldreich 1984). For critical cosmic density, fl = 1, it 
is typically comparable to the current Hubble time and the 
halo's spin axis is likely to slew by > 7° in 0.9 Gyr. 

A growing body of evidence suggests that fi is signifi- 
cantly less than unity, and one might be concerned that infall 
has largely dried up. In reality, it is not the global value of 
Q that is relevant here, but the effective value within the 
Local Group. Dynamical studies of the local group of ever 
greater sophistication all conclude that the mass of the Lo- 
cal Group is 4 - 8 x 10 12 M Q (e.g., Schmoldt & Saha 1998). 



If the halos of M31 and the Milky Way have flat rotation 
curves with amplitude v c , their masses interior to radius r 
are given by 

M _ / v c \ 2 / r \ 
10 12 M ~~ ' V220 kms-V V 100 kpc/ ' 

Consequently, the measured mass of the Local Group can be 
accommodated in the halos of these two dominant galaxies 
only if they touch each other. In this case, given that M31 
is approaching us, parts of M31's halo would now be falling 
into ours, and vice versa. Thus, even if all the mass of the 
Local Group had in the past been somehow organized into 
two quiescent halos, there would still be dynamically impor- 
tant infall at the present epoch. In reality, there is plenty of 
observational evidence that the Galactic halo is profoundly 
affected by infall even at r = 20 — 60kpc: the Galaxy is 
currently accreting both the Magellanic Clouds (e.g., Murai 
& Fujimoto 1980) and the Sagittarius dwarf galaxy (Ibata 
et al. 1997). Since observations of dwarf galaxies imply that 
both of these systems must be tracers of much more mas- 
sive streams of dark matter, the natural conclusion to draw 
is that most of the mass of the Local Group still lies outside 
the virialized halos of M31 and the Galaxy, and is falling in 
at a dynamically significant rate. 

This being so, it seems that warps probably are caused 
by infall-driven reorientation of the outer parts of virialized 
halos. 
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